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Abstract 

We consider a one-parameter family of short rate models which encompasses both Hull-White (normal) 
and Black-Karasinski (lognormal) models. We deduce a general form for the relevant Green’s function 
as an asymptotic; series, assuming only that the deviations of the short rate from the forward curve are 
on average small in absolute terms, and show how this solution can be parametrised in such a way as to 
fit the model to a term structure of zero coupon bond prices. We use the derived Green’s function to 
calculate conditional bond prices and pricing formulae for caps and floors to second order accuracy. The 
results are seen to take a form which is straightforward to compute using quadrature and even the first 
order expressions achieve highly favourable comparison with benchmark Monte Carlo computations for 
a wide range of market conditions with both long and short cap/floor maturities. 


1 Introduction 

Over a quarter of a century has elapsed since Black and Karasinski (1991) first proposed their eponymous 
lognormal short rate model for interest rates as an alternative to the normal model of Hull and White 
(1990) which, although relatively tractable, does not guarantee positive rates. As observed by Brigo and 
Mercurio (2006), the Black-Karasinski model turned out to be rather less tractable, which “renders the 
model calibration to market data more burdensome than in the Hull and White (1990) Gaussian model, 
since no analytical formulas for bonds are available.” More recently, Turfus (2016a) demonstrated how the 
distribution of conditional prices of zero coupon bonds under the Black-Karasinski model can be calculated 
consistent with a condition of no arbitrage, with errors which are third order in the short rate, with none 
of the assumptions about the smallness of volatility which are customarily made. Turfus (2016b) has shown 
how analytic expressions for caplet prices can also be derived for a family of models which includes both 
Hull-White and Black-Karasinski under an assumption of low volatility levels. This latter calculation was 
based on an asymptotic modelling approach to stochastic rates first propounded by Kim and Kunitomo 
(1999) in the context of vanilla equity option pricing. Their approach allowed them to address affine short 
rate models, in particular those of Hull and White (1990) and of Cox, Ingersoll and Ross (1985). 

In the present paper, we build on the above-mentioned results. We adopt the approach first taken by 
Hagan et al. (2005) in expounding their celebrated SABR model, namely to derive an asymptotically valid 
Green’s function for the underlying pricing PDE and to use this to deduce power series representations of 
option prices. Surprisingly little direct use has been made of this approach since then, although recently it 
has been revived and used by Pagliarani et al. (2011) in the context of local-stochastic volatility modelling. 
They have together with a number of collaborators applied their method inter alia to pricing Asian options 
on equity or other lognormal underlyings under a local volatility assumption (Foschi et ah, 2013). We follow 
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a similar approach in a short rate model context in deriving a second-order accurate asymptotic Green’s 
function for the Black-Karasinski and related models, then showing how this can be applied to bond and 
option pricing. 

We introduce in §2 a one-parameter family of beta blend models which includes Black-Karasinski and 
also the Hull-White model as special cases. In §3 we go on to derive a representation of the requisite 
Green’s function which is asymptotically valid for small deviations of the short rate from the forward curve, 
illustrating in the process how to calibrate the model to a term structure of zero coupon bond prices. 
The second-order accurate Green’s function obtained is then applied in §4 to the problem of calculating 
conditional bond prices to second order accuracy and further to obtain pricing formulae for caps/floors and 
capped/floored Libor payments. Favourable comparisons are found with Monte Carlo simulations in §5. 
Conclusions are stated in §6 and mention made of a number of ways the asymptotic approach proposed can 
be extended into a multi-factor modelling context to price hybrid derivatives. 


2 Modelling Assumptions 


We consider a beta blend short rate process depending on a parameter for the purpose of pricing interest rate 
derivatives. The parameter (3 may be taken to be constant or more generally time-dependent. This model was 
first proposed in the literature by Turfus (2016c) to describe the credit default intensity process in a hybrid 
equity-credit model and subsequently used also by Turfus (2016b) in an interest rate context. It encompasses 
both Hull and White (1990) and Black and Karasinski (1991) models as special cases, respectively for /? f 1 
and (3 = 0. 

In quantifying the short rate, we shall find it convenient to work with an ancillary process x t satisfying 
the following canonical Ornstein-Uhlenbeck processes: 

dx t = —a r (t)x t dt + u r (t) dW t (1) 

with dWt a Wiener process. Here the mean reversion parameter a r {-) : R + —> R + and the volatility 
oy(-) : M + —> M + are taken to be piecewise continuous functions. We follow the convention here and below 
that processes are distinguished from functions of t by use of a subscript t for the former. This ancillary 
process is related to the interest short rate r t by 


(1 ~/3(t))r t +(3(t)r(t) = (r{t) + (1 - /3(t))r*(t))£(x t , t) 


( 2 ) 


where we have defined £ (x, t) = £ (x, t, t) with 1 

£(x, t, u) := exp (Fp(u)<j> r (t, u)x — \Fp{u)4\.(t, u)I r ( 0, t)) , t> 0 
Fa(f \ , = (!-£(*)) f > n 

< j) r (t,u ) := e~ $<■ CXr ( s ') ds , u > t, 

I r (s,t) := f (J)l(u,t)(7r{u)du, t>s> 0 . 

J S 

Here r(-) : R + —» M + is assumed to be continuous, /?(•) : R + —> [0,1) is a parameter 2 controlling the skew 
of the model and 6 > 0 is a small parameter which can be introduced if needed to ensure positivity of the 
denominator. We further assume that xo = 0, with t = 0 the ‘as of’ date for which the model is calibrated. 
The function r*(t) is determined by calibration but will tend in the zero volatility limit to zero. The formal 
no-arbitrage constraint which determines this function is as follows 


(3) 

(4) 

(5) 

( 6 ) 


( 7 ) 


x It may be observed that £ (xt,t,u) thus defined is the stochastic exponential of Fp{t)cj> r (t,u)xt , t > 0. 
2 We shall see below how results for fH = 1 can be conveniently recovered by taking the limit as /3 —> 1~. 


E 


- f* ■ 

9 JO 


= D(0,t) 
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under the equivalent martingale measure for 0 < t < T m , where T m is the longest maturity date for which 
the model is calibrated, and 

D(t u t 2 ):=e-^ T ^ ds ( 8 ) 

is the ti-forward price of the ^-maturity zero coupon bond. 

We consider the (stochastic) time-f price of a European-style security which pays a cash amount P(xt) 
at maturity T, denoting this by f t = f(x t ,t). We note in particular that the price of a T-maturity zero 
coupon bond is obtained by taking P(x) = l. 3 We infer by standard means that in the general case f(x,t) 
satisfies the following backward diffusion equation: 4 

df ,, <9/1 2 , .d 2 f . . 

Si - aAt)x ai + r At) w= r{x ' t)f - (9) 


subject to the final condition lim f(x,t) = P{x), with r(x,t) := rt\ Xt =x and re-express this in operator 

t-yT~ 

notation as 

(J^ +£ ~ h{x,t)j f{x,t) = 0 (10) 

where 

1 Pi 1 f)f 

( 11 ) 

and h(x,t) := h(x,t,t) with 


h(x, t, u ) 


(r(u) + (1 — (3(u))r*(u))£(x, t, u) — r(u) 
1 - P{u) 


u> t, 


( 12 ) 


We further propose that the impact of h(.), considered as an operator, on /(.) is 0(e) in comparison to that 
of C and consider the limit as e —>• 0. Effectively we assume that the deviations of r t from the forward curve 
(or equivalently those of Libor rates) are small in absolute terms. 

Following Turfus (2016a), we note that satisfaction of the calibration condition (7) imposes conditions 
on r*(-) at 0(e 2 ) and higher. On that basis we propose an asymptotic expansion 


OO 


i=2 


(13) 


with r*{t) : 0{ei). Substituting for r*(-) into (12), we infer a similar expansion for h(x,t): 


h(x,t) = ^ hj(x,t), 

3 =1 


(14) 


with hj(x,t ) = O(e^). In particular 


hi(x,t) 


hj(x,t ) 


r{t)(£(x,t) - 1) 

1 - /3(f) 

r*(t)£(x,t), j> 2. 


(15) 

(16) 


In the following section, we will look to derive a Green’s function solution for (10) as an asymptotic power 
series in e. 

3 We will look to derive the functional form of /(•) implied by our model in this case and in the process to determine the 
conditions on r*(t) necessary to satisfy (7). 

4 Here the coefficient of / is simply the short rate rt conditional on xt = x. 
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3 Green’s Function 


3.1 Derivation of Green’s Function 

We will be interested in the so-called free-boundary Green’s function solution to (10) which tends to zero as 
x ±oo. We consider first the limiting problem with e = 0. The required limit solution can be written 


n , , , s n ,. , d (£-<j> r (t,v)x 

Go{x,t;£,v) = D(t,v)—N 


d£ V \/Ir{t,v) J' 


(17) 


where N(-) is a unit Gaussian cumulative probability distribution function. Guided by Kato (1995, see in 
particular Chapter IX, §2), we deduce the following result: 


Theorem 3.1 The function G(x,t;£,v) defined by the asymptotic series 

OO 

G{x, t; £, v) = ^2 G n (x, t; £, v) 

n— 0 

with Go(x,t; £,v) given by (17) and G n (x, t; £, v) = 0(e n ) defined recursively by 

n pV pO O 

G n (x,t;£,v) ^ / / G 0 {x,t;x 1 ,t 1 )h i (xi,t 1 )G n - i (x 1 ,t 1 ;^,v) dx^G, n> 1 

l = \ J t J — OO 

is a Greens function solution for (10) satisfying 

+ £ - h(x, t)j G(x, t\ £, v ) = 5(x - £)6(t - v) 

subject to G(x, t; t;, v) — > 0 as x ±oo for 0 < t < v. 

Proof 3.1 For the proof of Theorem 3.1, see Appendix A 
In particular for n = 1 we obtain from (19): 

n oo 

G 0 (x, t; x 1 ,t 1 )h 1 (x 1 ,t 1 )G 0 (x 1 ,t 1 ;tj, v ) dxidh 

-OO 

(£(x,t,ti)M tl - 1) G 0 (x,t;£,v) dt i, 


ft J —OO 

r(t i) 


Jt 1 — P(ti) 

where we have defined the shift operator 

M u G 0 (x, i; v) := G 0 (x, t; £ - Fp(u)(j) r (u, v)I r (t, u), v) 


(18) 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 
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and for n = 2: 


n oo 

G 0 (x, t; x 2 , t 2 )h 1 (x 2 , t 2 )Gi(x 2 , t 2 ; £, v) dx 2 dt 2 

-OO 

rv r oo 

G 0 (x, t; xi,ti)h 2 (xi,ti)G 0 (xi,ti-,£, v ) ofoqdfi 

(£(x,t,t 2 )M t2 - 1) [ —^f—(£(x,t,h)M tl -l)Go(x,t-,^,v)dhdt 2 

1 - P(t2) J t2 1 - P(*i) 


J t J —OO 

^ r(i 2 ) 


S (*2) 


1-^2) 


£(x,t,t 2 ) \ r ^) . £(x,t,h) _ i") 

1 — p(ti) V / 

M tl M t2 G 0 (x, t; £, v ) dhdt 2 


2 


^ J r 2 (h)£(x,t,h)M tl G 0 (x,t;£,v) dh 

(£(x,t,h)M tl - 1) (£(x,t,t 2 )M t2 - l)G 0 (a:,t;£,t>)cfticft 2 


^ r(i 2 ) 


r<2 F(ii) 


i-0(t 2 )J t i-0(h) 
r r(t 2 ) F( . , , r t2 r(t 1 ) 

-£(x,t,t 2 ) 


1 - /5(t 2 ) 


i -^£(x,t,h) ^e' F ' 3(tl) ' F ' 3( * 2 ^ r ' (tl ’ t2)/T ' (t,tl) - l) 

M tl Mt 2 G 0 (x, t; £, v) dhdt 2 


- J r 2 (h)£(x,t,h)M tl Go (x,t;£,v)dh, 


(23) 


where in the last expression we have changed the order of integration and reversed the roles of t\ and t 2 . 


3.2 Model Calibration 


Completion of the formal derivation of our asymptotically valid Green’s function requires that we ensure it 
is consistent with the 110 -arbitrage constraint (7), by choosing a suitable specification for the function r|(t). 
To this end we let fj be the time-f value of a unit cash flow paid at time T, whence the required final 
condition is frf = 1. Applying our second order accurate Green’s function for t = 0, we obtain 


D(0,T) 



( r(t 2 ) f t2 r(h) 

\l-0(t 2 )J o l-0(h) 


(tl)Fp (t2)<Pr(ti ,t2)I r {0,t!) 



We see that satisfaction of the required 110 -arbitrage condition = D(0,T) at O(e) is automatic on the 
basis of 110 0(e) term in r*(t) correctly having been assumed, whereas at 0(e 2 ) we must choose 


m 


r M f 1 (^1 { Ff,(u)Ff,{t)<t> r (u,t)Ir (0,u) _ J 

1 - 0(t) A 1 - 0(u) V J 


(24) 


We see that defining our model by the second order accurate Green’s function and choosing r*(t) = r%(t) 
will give us exact calibration. Note that in the limit as 0 f 1 we obtain in this way 

r*(t)—¥ f <t> r (u,t)I r (0,u)du (25) 

Jo 


as the well-known Hull-White limit. 

Finally, making use of (60) with (17), (21), (23) and (24), and supposing we can in general write 

OO 

p (x) =J2 p j(x), (26) 

j =0 
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( 27 ) 


with Pj{x) = 0(ei), we can compute a second order accurate representation of f(x,t) as 

2 2 —i «oo 

/ m~EE 

i =0 i=0 •' _ °° 

We note in passing that, under the assumed continuity of r(-), the form of our Green’s function imposes no 
smoothness restrictions on the target payoff function P(-), other than that it be integrable over any finite 
subset of R. 


4 Applications 


4.1 Zero Coupon Bond Pricing 

We consider next the conditional price /f of a T-maturity zero coupon bond using (27). In this case the 
payoff is given by P(x) = 1, whence we obtain 


/f = D(t, T) (1 - F f (x t , t) + F 2 t (x t , t) + 0(e 3 )) . 
where the F f(•) are 0{e l ) and given by 5 
rT r(fi) 


(28) 


(29) 


F i{x,t)=J t 1 (g(s,Mi) - 1) dtx 
F 2 (x,t) = (x,t) 2 + J £(x,t,t 2 ) 

1 -% 2 ) C (e F ' 3(tl)F/3(t2) ^ (tl ’ t2)/ ’' (Ml) - l) dh -r5(t 2 ))dt 2 . (30) 


In the Hull-White limit (/3 f 1) this becomes 

if ~ D(t,T) (1 - B*(t,T)(x + r 2 (f)) + h(x 2 - I r (0,t))B* 2 (t,T)) 

where 

B*(t,v) := J (j> r (t,u)du, 

which is consistent to second order with the exact result 

if = D(t, T) exp (—B*(t, T)(x + r 2 (f)) — C(0, t))B* 2 (t, T)) . 


(31) 

(32) 

(33) 


4.2 Caplet Pricing 

Consider a caplet which pays at time T a coupon based on the positive difference between tenor-r Libor and 
a strike K on a unit notional, based on a payment period [T — r, T\, which rate we shall denote L(t,T). For 
simplicity we assume no spread between forward Libor rates and the equivalent risk-free rates inferred from 
(2) above. However, it is not difficult to introduce an assumed deterministic spread by adjusting the value 
of the strike accordingly in the formulae derived below. We can express the caplet payoff at time T by 

Payoff T := max {(L(r, T) — K)S(T — t,T), 0} 

where f 2 ) is the day count faction calculated according to the relevant convention (typically actual/360 
or actual/365). We note in particular that, under our assumptions, the realised Libor rate is related to the 
stochastic zero coupon bond price /f calculated above by 

1 + L(t,T)5(T — t,T) = (frf-r)- 1 (34) 

5 We have here used the identity that, by symmetry, jf f(v) // f(u ) dudv = 1 ^jf f(u) du ^ . 
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whence 


Payoff T = max{(/f_ T ) 1 — (1 + KS(T — r, T)), 0} 

If we consider an equivalent-valued payment made at time T — r, this must be discounted by precisely the 
T-maturity zero coupon bond price observed at time T — r, whence we can write 

Payoff T _ T = kT 1 max {k — /^_ T , 0} (35) 


with 


1 


i + ks(t-t,t))‘ ^ 36) 

In other words we should consider a put option on the bond price. Writing the time-f price of this option as 
C T,K (x t ,t), we see this will satisfy (10) with final condition that 


C T ' K (x, T — t) = Payoff T _ 


t\X T -t—X- 


(37) 


We shall seek here to calculate the caplet price up to second order. The equivalent-valued payoff at time 
T — t is, using (28): 

P{x) = max {1 - k~ 1 D(T — r, T) (l - F?(x, T - r) + F?(x, T - r)) , 0} + 0(e 3 ). (38) 

We see that, since Ff (•) = 0(e), for the option to have non trivial value in the limit as e —> 0, we must have 
that 1 — k~ 1 D{T — t, T) = 0(e), so we make this assumption. Let us denote £* := inf{x | P(x) > 0}, in 
terms of which the Pj(x) in (26) can be written 


Po(x) = 0 

Pi (a-) = (1 - k~ 1 D{T - r, T)(l - Ff {x, T - r)))l x>e 
P 2 {x) = -n~ l D(T - t,T)F?(x,T- t) 1 x>( . 


Posing 


C T ’ K (x,t) = ^2CJ’ K (x,t) + 0(e 3 ) 
l=i 

in obvious notation and making use of (27), we obtain: 


C[' K 


(x,t) 


/ OO 

Go(x,t\£,T — r)Pi(^) df; 

-OO 


(D(t, T - t) - k-'DH, T))N(—di(£* — x,t,T — r)) 

+ K~ 1 D(t,T) f T 

Jt-t 1 - PW 

(£{x, t, u)N(—d 2 (£* — x,t,u,T — t)) — N(—di(£* — x,t,T — r))) du 


(39) 

(40) 

(41) 


(42) 


(43) 


where we define, for u, w < v, 


di(£,t,w) 


d 2 (£,t,u,w) 

d 2 {€,t,u,v,w) 

d 3 (^,t,u,v,w) 


$ 

VI r {t,w )’ 

di(£ — Fp(u)<j> r (u A w, u V w)I r [t, u Aw), t, w ), 
di(£ — Fp(v)<j) r (w, v)I r (u A w, w),t, w), 
d 2 (£ — Fp(v)(j) r (w, v)I r (u A w, w),t, u , w), 


(44) 

(45) 

(46) 

(47) 
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with the binary operators A and V denoting min and max, respectively. Likewise at second order, we have 

/ OO 

(Gr Or, P,£,T- r)Pr(0 + G 0 (x, i; T - t)P 2 (0K 

-OO 

r t ~ t y(u) f°° 

= -(1 - k~ 1 D(T-t,T)) / --— (£(x,t,u)M u - 1) / G 0 (x,t;£,T - t) d£du 

Jt 1 - p(f) df* 


It 1 - /3(f) 
"(f) 

-/3(f) 


r r(u) f°° 

-k~ 1 D(T-t,T) -- ——(£(x,t,u)M u - 1) / G 0 (x,t;£,T-T) 

Jt 1 - PW d ? * 

(£, T — t) d£du 

nOO 

— n~ 1 D(t, T) J * Go(x, t\£,T — t)F 2 (£, T ~ r)dt; 

f T-r 

= -(D(t,T-T)-K~ 1 D(t,T)) j - 


r{u) 

-/3(f) 


(£(x, t, w)iV(—d 2 (£* — x,t,u,T — r)) — N(—di (£* — x,t,T — r))) du 


dr-r 1 - /3(f) 


t 1 - /3(f) 

^/j(i*)^(®)«r(»A(T-r),«)J P (t,uA(r-T))£( a . } t} u )N(-d 3 (C _ ajj t> U,V,T — r)) 

— _/V(—d/|(£* ~ ti u,v,T — r))^ dudv 


l D(t,T) 


"(f) 


T—r 1 ^ /3(f) 


£{x, t, v) 


r r(u) 


1 - £(«) 


(£(x, t, w)./V(—d 2 (£* — x,t,u,T — r)) — N(—di(£* — x, t,T — r))) dudv 

+ n~ 1 D(t,T) f r 2 (v)£(x,t,v)N(—d 2 (f;* — x,t,v,T — T ))dv. (48) 

Jt-t 


Setting x = t = 0 and combining terms, we obtain: 

G T ’ K (0,0) - (£>(0, T - r) - k _ 1 D(0, T)) JV(-di (f, 0, T - r)) 


+ K - i D(0,T) 


r{u) 


Jt—t 1 - /3(f) 

-(D(0,T-t)-k- 1 £>(0,T)) [ 

Jo 


(N(-d 2 (CA U, T-t))- iV(-di(r, 0, T - r))) du 


r(u) 


o 1 - /3(f) 

(ZV(-d 2 (f,0, u,T — r)) - /V(-dr(r, 0, T-t))) du 


- k^D^T) 


r(v) 


r(u) 


0 F /3 (u)Fp (y)<f> r (uA(T—r) ,t>) J r (0,uA (T—t)) 


' T _r 1-/3(^) Vo l-/3(w) 

(N(-d 3 (C,0,u,v,T - t)) - N(—d 2 (£*,0,v,T — r))) dudu 


+ K - 1 D(0,T) [ T f ^|L(7V(-d; ( r,0,f,f,T-r)) 

Jt-t 1 - /3(f) do 1 - /3(f) 


/t-t 1 - /3(f) do 1 - /3(f) 

- iV(-d 2 (r, 0, v,T — T)) + N(-d 2 (C, o, u,T — T)) - /V(-di(r, 0, T - r))) dudv, (49) 


with G(e 3 ) errors, where the first two lines constitute the first order solution and the remainder provides the 
necessary second order adjustment. A resemblance to the standard Black formula is here evident in the first 
order solution, the most notable difference being that the d 2 (-) argument is sampled across the integration 
range u £ [T — t,T ], rather than being considered fixed, as with the corresponding (exact) Hull-White result. 

We see that the integrand in (49), although it is not defined for 3 = 1, is yet well-behaved in the limit as 
/? f 1, which it is a straightforward matter to compute. However, given the existence of an exact solution, 
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this is of less interest than when /3 < 1. The Blaek-Karasinski result is obtained by setting /3 = 0 and Fg = 1 
in (49); also in (29) and (30), which are needed to obtain £* from (38). 


4.3 Floorlet Pricing 

By extension, the corresponding second-order accurate expression for the floorlet is 
F TJ< (0,0) - (« -1 -D(0, T) — D(0,T — r)) N(d, (£*, 0, T - r)) 

-k-'D^T) [ 1 r( ^ . (N(d 2 (C,0,u,T-T))-N(d 1 (C,0,T-T)))du 


IT-t 1 — ft( u ) 


T-r - 


— (k _1 D(0, T) — D(0,T ~ t)) - 


(u) 


-P(u) 


+ k~ 1 D(0,T) 


- K~ i D(0,T) 


*(v) 


r r(u) 


(• N(d 2 (f ,0,u,T — t)) — N(di(£*, 0, T — t))) du 

Fp(u)Fp (v)(f> r (uA(T—T),v)I r (0,uA(T—r)) 


T—r 1 - Vo 1 - £(u) 

{N(d 3 (£*,0,u,v,T - t)) - N(d 2 (C,0,v,T - r)))dudv 


r(v) 


r(u ) 


T—r 1 - £(*0 do 1 - /?(«) 


(Jv(^(r,o,«,v,r-r)) 


- N(d 2 (C, 0, v,T — r)) + N(d 2 (C, 0, t», T-r)) - W(di(f*, 0, T - r))) dwdu, (50) 

with the definition of £* unchanged (although technically the payoff function now has the opposite sign, 
whence £* will be a supremum for x rather than an infimum). It is also readily verified that use of (49) and 
(50) leads to 


C T ’ K { 0,0) - F t ’ k ( 0,0) = D{ 0, T) r{u) du - KS{T - r, T)j , 

so put-call parity is seen to be exactly satisfied under our second order approximation. 


(51) 


4.4 Capped/Floored Libor Flows 

We consider also briefly the pricing of Libor flows which are capped and/or floored. We note that the impact 
of a cap at level Kq on a Libor flow is by replication given by subtracting the value of a caplet with strike 
Kc ■ Since the time-t value of a Libor flow struck at T — t > t for payment at T is, using (28), 

L T (x t ,t) = fr T -f?, (52) 

the capped value will be 

CL t,Kc (x t ,t) = L T (x t ,t) - C T ’ Kc (x t ,t). (53) 

and the corresponding value of a Libor flow floored at Kp will be 

FL T ’ Kf (x t , t) = L T (x t , t) + C T ’ Kf (x t , t). (54) 


If both a cap and a floor are applied, the result will be 

CFL T ’ Kc ’ KF (x tl t ) = CL T ’ Kc (x t ,t) + C t ' Kf (x t ,t). (55) 


5 Comparison of Results 

In the following we focus on the pricing of caps: given the demonstration of put-call parity in (51), the error 
associated with the asymptotic approximation in the corresponding floors is guaranteed to be identical, so 
does not merit separate investigation. 
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Figure 1: Black-Karasinski prices for 5 y maturity cap with 6m Libor tenor 
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Figure 2: Blend prices with /? = 0.5 for 5 y maturity cap with 6m Libor tenor 


Calculations were mainly done for the important Black-Karasinski case, using a first order asymptotic 
representation. That is, we ignored the contribution from in (38) and used only the first two lines of (49). 
Comparisons were made with a Monte Carlo simulation using a conservative 500,000 simulations, ensuring 
any Monte Carlo error did not impact significantly on the quality of the comparison. The short rate is 
taken to follow the profile of the USD market for January 2017 rising from around 1% to around 3% over 
a five year period. The (lognormal) local volatility is taken to be 30%, with a (constant) mean reversion 
rate of 25%. The results are illustrated in Fig. 1 for various moneyness levels. As can be seen, they are 
indistinguishable to visual accuracy. It was observed that in-the-money caplets are slightly underpriced and 
the out-of-the-money caplets correspondingly overpriced, but the latter effect is considerably more marked, 
resulting in almost all caps being to some degree overpriced. 

By way of comparison, we also consider a case with /3 = 0.5. On this occasion we chose a local volatility 
of 4%, which yields comparable prices. Results are illustrated in Fig. 2. The level of error is seen in general 
to be less than that for the Black-Karasinski case, with underpricing again found for in-the-money caplets 
and overpricing for out-of-the money caplets, the latter effect being the stronger. 

Since the Black-Karasinski model is the case of greatest interest, we focus in the remaining tests on 
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Figure 3: Black-Karasinski prices for caps with different maturities 
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Figure 4: Dependence of Black-Karasinski prices for 5y caps on volatility level 


that case. We consider next how the asymptotic, expansion performs in the case of different maturities up 
to 20y. For the purposes of allowing comparative analysis across different maturities, we focus on ATM 
options by moving the strike of the caps appropriately. Results are shown in Fig. 3. They remain fairly 
indistinguishable to visual accuracy. As is evident, the error increases with maturity but of course so does 
the cap value. Overall the relative error increases, so the approximation does deteriorate slightly for longer 
maturities. 

Next we consider the variation of cap prices with the volatility a r (t). Results are shown in Fig. 4. As 
with time to maturity, the error is seen to grow as volatility increases, as does the price, both approximately 
linearly in the high volatility limit. Even at the very high level of 50% volatility the error is still only about 
0.4% in relative terms. 

A similar comparison is illustrated in Fig. 5 of the dependency of 5y cap prices on the assumed volatility 
mean reversion rate a r (t). Not surprisingly, as the mean reversion rate approaches zero allowing the volatility 
to grow without bound as T increases, the errors in the asymptotic approximation become larger, but even 
in such circumstances no greater in relative terms than those observed in Fig. 4 above. 

We consider the impact of the Libor tenor on 5y cap prices in Fig. 6. As is to be expected, the tenor does 
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Figure 6: Dependence of Black-Karasinski prices for 5y caps on Libor tenor 


not impact greatly on the price. But it does impact on the error, with the caps based on ly Libor resulting 
in significantly smaller discrepancies. With the magnified y-axis in this case, the slight overpricing of the 
asymptotic model is visually detectable. 

Given that our asymptotic expansion is based on an assumption of small fluctuations in the rates in 
absolute terms, it is to be expected that the accuracy will be lessened as the level of the yield curve is raised 
with fixed lognormal volatility levels. The impact of such a change is considered in Fig. 7, where the yield 
curve is parallel-bumped upward by up to 8%. To enable a fair comparison, the strike is raised by the same 
amount in each case, ensuring that options remained roughly at the money. As can be seen, the option value 
increases due to the (fixed) lognormal volatility having greater PV impact for higher rates. As expected, 
the error also increases, roughly in proportion. Even in the fairly extreme case of a 10% interest rate, the 
relative error in the 5y cap PV was seen to be only 0.5% in relative terms. In circumstances where rates 
exceeded the 10% level, it might be considered worthwhile to include second order terms in the calculation 
so as to ensure greater accuracy. 

In conclusion, we infer that even our rather simple first order approximations ought to be usable for 
many practical purposes, particularly calibration and risk management where exact reproduction of market 
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Figure 7: Dependence of ATM Black-Karasinski prices for 5y caps on IR curve level 


prices is less essential than in pricing for trading purposes. If greater accuracy is required, particularly in 
the Black-Karasinski case, the inclusion of second order terms ought easily to provide enough accuracy with 
very little additional computational burden. 

6 Conclusions 

A new family of short rate models is introduced which offers a blend between the well-known Hull-White 
and Black-Karasinski models. An approximate Green’s function was derived in §3 for the resultant pricing 
PDE, based on the assumption of the deviations of the Libor rates from their forward values being small in 
absolute terms. This Green’s function was used to derive the second-order accurate expression (28) for the 
conditional price of zero coupon bonds. Further, second-order accurate expressions were derived for caplet 
and floorlet prices, which results can be combined to allow analytic pricing of caps and floors, or indeed 
of capped/floored Libor flows. The caplet price is given by (49) and the corresponding floorlet price by 
(50). These expressions are all straightforward to compute. Similar results can be found for the Hull-White 
case but are of less interest on account of availability of an exact analytic solution. Even the first order 
expressions were found to yield highly favourable comparison with benchmark Monte Carlo computations 
for a wide range of market conditions with both long and short cap/floor maturities. 

The expressions for (3 < 1 should prove useful in calibrating the model to at-the-money market prices for 
caps and floors, which task would otherwise be computationally intensive. However an additional important 
use of the approach presented here is in the context of the short rate model being combined with another 
component model (credit intensity, equity, FX and/or another interest rate, say) for hybrid derivative 
pricing purposes. Of particular interest is use of the Black-Karasinski version of the model to describe 
credit intensities since this guarantees their positivity. Such work is currently under way and numerous 
preliminary results have been obtained. See for example Turfus (2017a) and Turfus (2018b) for application 
to the valuation of credit default swaps taking account of wrong-way correlation risk, Turfus (2017b) for the 
pricing of capped/floored Libor payments in a quanto CDS pricing context, Turfus (2017c) for the pricing 
of CDS on an underlying interest rate swap, Turfus (2017d) for an application to calculating the CVA on 
a collateralised portfolio of equity options, Turfus (2017e) for the pricing of convertible bonds and Turfus 
(2018a) for an example involving the pricing of quanto CDS and CDS with a quanto loss cap. 

Much work has also been done by other authors using perturbation approaches to obtain analytic 
approximations for prices of various option types under local and/or stochastic volatility modelling 
assumptions. See, for example, Pagliarani et al. (2011), who considered equity option pricing under a 
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local volatility assumption, obtaining a perturbation expansion for the relevant Green’s function much as 
we did here, and using it to derive asymptotic expressions for option prices. Their approach was applied 
also to Asian option pricing in Foschi et al. (2013) and extended, with the use of some Fourier analysis, to 
incorporate Levy jumps in the dynamics of the spot underlying in Pagliarani and Pascucci (2013). A review 
of a number of other papers which have presented asymptotic option pricing formulae in recent years has 
been given by Turfus and Schubert (2017). An interesting prospect for future work would be to extend the 
above by applying the asymptotic short rate model for stochastic rates/credit described herein alongside 
a local-stochastic volatility model to capture their combined impact, as was done by Funahashi (2015) for 
European equity options. 


A Proof of Theorem 3.1 

We demonstrate that (18) is a solution to (20). We start by defining 

N 

G N {x,t-,£,v) = ^2 G n (x,t;£,v). (56) 

n—0 

Using (19) and the fact that by definition 

(J^ + £ ) G °^’ t; ^ V ^ = S ^ X ~ ^ ( 57 ) 

we see that 

d n 

—G n (x, t; g, v) = hj(x, t)G n -j(x, t; g, v) - CG n (x, t; g, v), n > 1, (58) 

i=i 

whence, using (56) and (57), 

( a \ N n 

— +CjG N (x, t; £, v ) = S(x - £)5(t - v) + EE hi(x, t)G n -i(x, t; £, v) 

' n—1 i =1 

N N-i 

= 5>(*,t) E G r {x,t;£,v) + 5(x - £)5(t - v). (59) 

i—1 r—0 

Finally, taking G(x,t-,f;,v) to be limjv^oo G N (x, t; £, v) and using (56), (59) and (14) we obtain 

idt + G ^ X ' t; ^ V ^ = v ) + S ( x ~ “ u )> ( 60 ) 

confirming that G(-) so defined is the required Green’s function. This concludes the proof of the theorem. 
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